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A crystal  of  CVD  zinc  selenide  was  placed  in  the  steady  2 kW  beam  of  a 
10.6  fjjn  laser.  The  beam  was  approximately  gaussian  with  l/e  radius  of 
0.  24  cm.  E’ifty  frame  per  second  motion  pictures  were  taken  of  the  time 
dependent  optical  distortion  as  observed  with  He-Ne  light  and  a Mach- 
Zehnder  interferometer.  Thirty-one  cir cular^ringes  originated  in  10 
seconds. 
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theoretical  solution  is  derived  for  the  heating  rate,  rate  of  development 
of  thermoelastic  stress  and  rate  of  change  of  optical  path  as  a function  of 
time,  radius,  and  surface  cooling  rate.  The  expressions  are  obtained  in 
closed  form  for  a gaussian  laser  beam  and  presented  in  non-dimensional 
form.  The  theory  is  then  compared  with  the  experiment.  For  the  first 
five  seconds  of  heating,  there  is  good  agreement.  In  the  next  four  seconds, 
the  heating  rate  is  higher  than  calculated  and  the  discussion  contains  some 
possible  explanations. 
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TEMPERATURE,  STRESS  AND  OPTICAL  PHASE  DIFFERENCE 
FOR  A WINDOW  TRANSMITTING  A CW  LASER  REAM 


The  transmission  of  a laser  beam  through  a plane  window  gives  rise 
to  a change  in  optical  path  length  as  a result  of  three  physical  phenomena; 

(1)  change  in  window  thickness  by  heating  and  thermal  expansion;  (2)  change 
in  index  of  refraction  as  a result  of  heating  and  (3)  change  in  both  thickness 
and  index  as  a result  of  thermoelastic  stress  set  up  by  temperature  gradients 
during  the  time  of  illumination  or  subsequent  cooling  of  the  window.  In  this 
paper,  a solution  to  the  optical  path  difference  (OPD)  resulting  from  these 
three  mechanisms  is  derived  in  analytic  form  and  the  results  applied  to  a 
window  experiment  where  the  time  dependent  OPD  for  a plane  window  was 
experimentally  determined  during  and  after  illumination  by  a laser  beam. 

The  problem  of  the  time  dependent  heating  of  a laser  v/indow  by  an 
isotropic  and  centrally  positioned  laser  beam  has  already  been  solved  by 
various  methods.  In  the  present  work  the  Fourier  Bessel  trans- 

form is  used  and  it  yields  the  temperature,  strain,  and  stress  distributions 
in  closed  form.  The  window  response  can  thus  be  evaluated  without  the 
need  for  numerical  integration  methods,  at  least  for  some  problems  of 
practical  importance. 

A typical  laser  window  material  is  isotropic  with  failure  characteri- 
stics generally  described  as  brittle.  For  an  acceptable  window,  the  absorption 
of  laser  radiation  in  it  is  sufficiently  low  that  the  transmitted  beam  may  be 
approximated  by  an  intensity  distribution  independent  of  z.  The  window 
has  a radius  b,  a thickness,  f.  The  incident  laser  beam  that  is  used  is 
normal  to  the  window  surface  and  has  a gaussian  intensity  distribution 
given  by 

? 2 2 

1 = V a e'^  ^ /n  (1) 


The  first  two  authors  use  finite  Fourier  Bessel  transforms  to  state  the 
problem  but  use  numerical  methods  for  the  integrations.  The  third  author 
does  not  employ  a gaussian  laser  beam. 
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where  P = laser  power  and  a is  a constant.  The  window  has  an  absorption 
coefficient  of  P per  cm  (Pf  « 1)  so  that  the  heating  rate  in  degrees  per 
cubic  cm  per  second  is  given  by 

2 2 

T = T e"^  (Z) 

o 

and 

2 


f3) 


o Trpc 

where  pc  is  the  specific  heat  per  unit  volume  of  the  laser  window  (assumed 

to  be  a constant).  Thus  far  nothing  has  been  said  regarding  reflections 

from  the  incident  and  emerging  surfaces  of  the  window.  This  should  be 

considered  because  P is  defined  in  terms  of  the  intensity  of  radiation 

inside  a medium.  If  one  calculates  the  reduction  in  intensity  resulting 

from  surface  reflections,  and  also  the  increase  in  absorption  resulting 

from  multiple  internal  reflections  of  the  transmitted  radiation,  the  two 

(4) 

effects  compensate  each  other.  One  finds  that  the  absorbed  intensity 
A « pf  independent  of  the  reflection  coefficient  and  the  reflected  intensity 
R « 2 independent  of  P (Rj  is  the  reflection  coefficient  from  a single 

optical  surface).  Thus,  nature  has  cooperated  in  this  instance  in  making 
it  unnecessary  to  formally  include  the  micro  details  of  the  reflection  and 
transmission  process.  All  that  is  required  in  Eq.  (3)  is  that  pf  be  small 
and  that  the  two  surfaces  of  the  window  be  not  accurately  perpendicular 
to  the  wave  front  in  order  that  the  localized  interference  fringes  resulting 
from  internal  reflections  are  spaced  much  closer  than  l/a  apart. 

The  equation  leading  to  the  temperature  distribution  in  the  window, 
as  described,  is  the  Fourier  Heat  Conduction  Law. 


2 

a^T 
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K dt  K 
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0 < r < b 
0 < z :S  £ 


where  k is  the  thermal  diffusivity  and  the  last  term  is  the  heat  generated 

per  unit  volume.  The  boundary  conditions  are  that  -5 — =0  at  z = 0 and  f 

{)  T ^ ^ 

and  = 0 at  r = 0 and  b.  The  case  where  surface  cooling  occurs  will 

be  treated  subsequently  as  an  extension  of  the  solution  of  Eq.  (4). 


The  differential  equation  is  solved  by  two  calculations  in  which  the 

temperature  at  the  center  is  first  a step  function  and  then  the  solution  is 

rewritten  for  a continuous  central  temperature.  The  first  step  makes 

use  of  the  Hankel  transform  of  order  zero  (also  called  the  Fourier  Bessel 

transform)  T (^)  = r J^fr^)  T (r,t)  dr. 

o 


T(e)  = A 


r 

° Jo 


r (r^)  exp  (-a^r^)  dr 


(5) 


where  is  the  temperature  at  the  center  of  the  window,  r = 0,  The  inte- 
gration is  given  in  closed  form  by^^^ 

f(4)  - (A^Za^)  exp  (-^^/4a^) 

After  integration,  the  solution  is  found  by  taking  the  inverse  transform 


T(r,.,  = J 

f e'*^  ^ T(^)  e (^r)  de 

o 

(6) 

X 

^-4^  [l  + 4 K a^t]  /4a^^  % (|  r) 

(7) 

T(r,t)  = + 4 K a^t 

]^exp^-r^a^/[l  + 4k  a^t]^ 

(8) 

The  second  calculation 

4.V  (6) 

theorem 

uses  the  Duhamel  superposition 

Ad 

and  since 

f A^  (t)  ~ T (r,  t - T)  dr 
o 

A (t)  = f 

T dT  = T fdT  = r (t)  . 

o / 

Jo 

o oj  o 

T(r,t)  = tC 

(7) 


t 


exp  I (-r^a^/[l+4Ka‘^ 


1+4  Ka  (t-  T) 


This  integrates  to 
T(r,  t)  = 


(10) 


(in 
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Here  Ei(x)  is  the  exponential  integral  of  negative  arguments,  tabulated 

/ gv 

in  Jahnke  and  Emde' 

(5) 


as  the  function  Ei(x)  = - f~  t ^ exp  f-t)  dt. 

. . /x 

(Abramowitz 

Equation  (11)  gives  in  closed  form,  the  time  and  radial  dependence 
of  the  heating  of  a window  by  a laser  beam  of  gaussian  intensity  profile 
when  the  window  is  larger  than  the  laser  beam,  there  is  negligible  heat 
loss  from  surfaces  and  edges  and  there  is  negligible  reduction  in  laser 
intensity  as  a result  of  internal  absorption.  Figure  1 is  a plot  of  the  non- 
dimensional  temperature  rise  T/T^  t against  the  non-dimensional  window 
radius  ra,  as  given  in  Eq.  (11). 

Window  Stress  and  Strain 

The  stress  and  strain  in  a window  transmitting  a laser  beam  of 

gaussian  intensity  can  also  be  stated  in  closed  form.  The  basic  equations 

(9) 

are  those  for  solid  cylinders.  If  the  window  is  thicker  than  the  laser 
beam  mean  radius,  one  should  use  the  long  cylinder  expressions  (plane 
strain)  rather  than  the  thin  plate  formulae  (plane  stress)  because  the  latter 
theory  ignores  axial  stress.  Since  failure  of  brittle  materials  occurs  where 
the  shear  is  a maximum,  the  thin  plate  theory  leads  to  the  conclusion  that 
windows  are  stronger  than  they,  in  fact,  are  because  the  axial  stress  terms 
are  the  largest  of  the  principle  stresses.  The  equations  for  plane  strain 
are; 
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( l-v  ) a /aE  = 


( 1-v  ) <j^Ae 


Trdr 
b 


(l/i 


■'I 


Trdr 


{\/\T)  f Trdr  + (l/r  )T  Trdr 

Jo  Jo 


(l-v)CT„/aE  = (Zv /'O')  J Trdr  - T 


(12) 


(13) 


(14) 


In  the  above,  a is  the  linear  coefficient  of  expansion,  E is  the 
elastic  modulus,  v is  Poisson' s ratio;  o is  strain  and  2b  is  the  diameter 
of  the  window.  These  expressions  contain  ^Trdr  which  fortunately  is 
integrable  from  Eq.  (11). 
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(151 


The Urdr  has  a much  simpler  physical  significance.  It  is 
proportional  to  the  total  energy  absorbed  in  the  window  and  therefore  can 
be  related  to  the  incident  laser  power. 

- - /1,2  ry.  (16) 


jT  'T^dr  = ^ = (b  /2)  Tj 


where  is  the  temperature  rise  for  infinitely  high  window  conductivity 
and  is  a constant  for  a particular  experiment. 

Upon  substituting  the  expressions  (15)  and  (1 6)  into  the  stress 
Eqs.  (12)  to  (14)  it  follows  that 
2( 


1:11%  v/_i B_\ 

a ET  - T \ 2.  2 ~ 2 2 f 

\ a b a r / 


2(  l-v)CT 


a ET 
2(  l-v  )a 


where 


T t 

aET  T 2^2  ' ^ 

a b 


exp  (-r^a^)  * ( i+4/ca^t)  exp  j- 


(17) 


(18) 


(19) 


B 


4Ka^t 


V l+4xa^t/ 


Figures  2,  3 and  4 are  plots  of  the  non-dimensional  radial,  tangential  and 
axial  stresses  respectively 


*r,  9,  z 
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VALUES  OF 
Afco^t 


T 


T 


From  which  ^h.e  individual  stress  valuos  may  be  obtained  for  an  experi- 
ment by 

S . fU’aE 
r,  0,  T. 

B,  z SiTKpc  ( I - V ) 

li  ran  be  seen  from  these  carves  that  the  arial  stress  is  about  twice  as 
l.ir^e  as  the  others. 

The  axial  strain  for  a thin  cylinder  (plane  stress)  ie  given  by 


C = a T + cr  /E  - v (o  + a*)/E 

Z 2 ^ 


The  plane  stress  equations  are 


= aE  j^Tj/2  - (l/r^) 

= qE  j^Tj/2  + (''■■)/  Trdr  * 


(22) 


a = 0. 
z 


(23) 


Tliese  lead  to  the  expression  for  axial  strain  when  the  laser  beam  nearly 

completely  fills  the  window 

e = oT  -i  va  (T  - T.) 
z 1 

For  cylinders  that  are  thicker  than  the  laser  beam  radius  one  uses  the 
stress  Eqs.  (12)  to  (14)  but  in  addition,  one  should  rewrite  Eq.  (22)  by 
replacing  EbyE/(l-v^);v  by  v /(1-v)  and  a by  a (1+v).^*^^  This  leads 
to  the  thick  window  expression  for  axial  strain. 

'z  = “ TT^  “ <T  - T,)]  -24) 

In  terms  of  non-dimensional  quantities,  one  has  the  non-dirr.ensional  axial 
strain 


E 


z 


4 7T ( 1 - V ) PCK 

ntv)2 


a b 


T 

T t 
o 


(4/ta  t) 


(25) 
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In  early  stages  of  window  heating,  the  second  term  in  the  bracket  may  be 
neglected  compared  to  1 so  that 


E,  ~ (4Ka^t)  f26) 

^ T t 
o 

This  is  plotted  in  Figure  5 for  various  values  of  the  non-dimensional  radius 
and  transmission  time. 

F ailure 

In  triaxial  loading  of  brittle  materials,  failure  occurs  when  the 
maximum  shear  stress  equals  half  the  yield  stress,  as  determined  in 

a uniaxial  test. 

There  are  three  maxima  of  shear  which  occur  at  45  degrees  to  the 
principle  stress  directions.  We  expect  failure  when 


a_)  = -r  a 
2 2 y 


where  and  are  any  of  the  three  principle  stresses  taken  in  whatever 
order  gives  the  largest  r,  Since  the  axial  stress  is  approximately  twice 
the  other  stresses,  failure  is  predicted  when  « 2o^  and  this  first  occurs 
on  the  central  axis  of  the  laser  beam.  The  max  shear  stress  is  plotted  in 
Figure  6. 

Optical  path  changes  in  a plane  window  that  is  heated  by  a laser 
beam  v/as  observed  by  including  the  window  in  one  leg  of  a Mach  Zehnder 

(4) 

interferometer.  The  dependence  of  the  OPD  on  the  temperature  is 
found  by  evaluating  the  individual  contributions  of  several  phenomena, 

T T 


OPD  = An£  = 


X " 


djg 
s dT 


dT  + 


dT 


<21) 


where  n^  is  the  index  of  refraction  as  customarily  measured  at  zero 

stress,  f is  the  window  thickness  and  C is  the  thickness  measured  at 

s 

zero  stress.  Through  the  use  of  Eq.  (24), 


dT  UT 


1-  V 


a 

s 


(28) 


-I  1- 


] 


Here,  is  the  coefficient  of  expansion  as  measured  at  zero  stress.  It  is 
to  be  noted  that  Eq.  (28)  already  includes  the  thermal  expansion  term  even 
in  the  absence  of  temperature  induced  strains.  Also, 


dn 

dT 


Q da 
o n r 

do  dT 


da 


dn  

3a„  dT 


B 


9 n 


da 


dn 


da  dT 


+ 


dT 


(29) 


B z 

The  last  term  in  Eq.  (29)  is  required  because  there  is  a change  in  index 
of  refraction  with  temperature  even  when  the  temperature  is  uniform 
and  there  is  zero  thermoelastic  stress.  Following  the  usual  classical 
derivation  of  stress  birefringence,  one  has  for  cubic  crystals^  trans- 
mitting a plane  wave  in  the  z direction, 


1 


2 


9 n 
0 a 


do 


( iO) 


He 


9 n 


re  q,,  and  q._  are  the  stress  optic  coefficients,  q,  , =— and 

^11  ^12  ^ 119a 


_ 9 n 

'l2  9a 


8 n 
9 a ,, 


dn 

dT 


7 t" 


da 

z 

dT 


d 

*^12  dT 


(a^  + or^) 


T dn 

r~ 


s 

dT 


(31) 


The  derivatives  of  the  stresses  are  obtained  from  Eqs.  (17)  to  (19). 
The  evaluation  of  these  stress-optic  terms  is  very  simple  for  cubic 
crystals  because  the  B term  is  not  present  in  either  a or  in  the  sum 

(»,  tag). 


da 

z 

dT 


d(a^  + Oq) 

dT 


uE 

1-v 


OPD  = 


Tjg 

1-v 


P 2a  E 

(l+v)'^Qg- ^ (qjj  + q^2^  + (1  - v)  ® 

n 

s 


dn 
£ 

dT 


(32) 


Note  that  the  last  term  in  the  bracket  of  Eq,  (32)  is  often  the 
dominant  term. 
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Surface  cooling  of  windows  through  the  use  of  a flushing  gas  stream 
is  an  effective  method  of  removing  laser  beam  generated  heat.  This  is 
particularly  useful  for  window's  w'ith  antireflection  coatings  because  a large 
part  of  the  window  heating  may  originate  in  the  surface  and  it  is  efficient 
to  remove  it  at  the  source. 

Equation  (4)  is  now  expanded  to  include  a term  for  removal  of  heat 

2 

from  two  surfaces  with  a surface  heat  transfer  coefficient,  h W/cm  deg. 


1 AZ 

r Or 


i ZZ 

K 9 1 


T ZhT 
k' 


The  solution  to  this  equation  proceeds  in  a similar  way  to  that  for  Eq,  (4); 
that  is,  by  first  finding  a solution  for  a steady  step  addition  of  heat  by  a 
gaussian  source,  (t/k  = 0 in  Eq.  (3  3))  and  then  applying  the  Duhamel 
superposition  integral  to  calculate  the  result  for  a time  dependent,  con- 


tinuous heat  source. 


Equation  (33)  becomes 


i ZZ 

r 8 r 


i.  ZZ 

K 9 1 


= h’T 


where  h'  = 2h/k£.  The  solution  for  the  step  heat  addition,  corresponding 
to  Eq,  (8)  is 

T^(r,t)  = I^T  /(I  + 4 Ka^t)J  exp  [^-/ch't]  exp  ^-r^a^/(l  + 4 j(a^t)j  13 

In  the  above  T is  a temperature  change  at  r = 0 which  will  be  evaluated 
shortly.  This  equation  is  now  inserted  in  Eq.  (9) 


T(r,  t)  = T 


L-Kh'  (t-  T)J 


2 2 

L-r  a ; 


^l  + /fa  (t-r) 


1+4  Ka  (t-  T) 


The  integration  is  performed  by  expanding  the  first  exponential  in  a power 

series  exp[-  Kh'  ( t - 7)]  ^ 1 - Kh'  (t-T)+.,..  For  the  expe  riment  that  we 

will  try  to  explain,  the  Kh' t •^<  1 for  t « 100  sec  and  therefore  only  the  first 

two  terms  of  the  expansion  are  sufficient  at  early  times.  This  gives  two 

terms  to  the  integral  of  Eq.  (36)  T(r,  t)  = + T^.  is  identical  to 

Eq.  (10)  for  which  the  solution  is  stated  in  Eq.  (11).  This  also  permits 

identification  of  'i'  as  identical  to  T . 

o 

T^  can  be  found  in  closed  form  by  integration  by  parts. 
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T^'r.t) 


T - I exp  F-r^a^l  - (1  + 4 Ka^t)  exp  f^-r^a^/(  1 + 4 fca^t) 

° (4/fa  I t-  J L 


2 2 

+ (1  + r a 


) Ei  -Ei  [- 


2 ^ /i  1 j_  A 
r a /( 1 + 4 Ka  t) 


) 


The  complete  solution  of  Eq,  (36)  is  thus  obtained  by  adding  Eq.  (11) 

t 


T(r,t)  = 


[• 


I + 


h'  (1+r^a^) 
4a^ 


Ei  (-r^a^)  - Ei 


\ I +4(ca^t  /J 

(-tSt)]) 


+ I exp  ( - r^a“)  - ( 1 + 4 /ca^t)  exp  I - 

4a 

where  «h'  t « 1.  On  the  axis,  (r  = 0)  this  equation  must  be  restated  to 
remove  the  indeterminate  form. 


(.37) 


Urn  T(r,  t)  = 
r=0  4ks 


^ In  ( 1 + 4 Ka^t)  ( 1 4 \ - 4Ka^t 

a 1_  V 4a  / . 


(38) 


Comparison  with  an  experiment.  In  an  experiment,  using  a CW  laser 

(4) 

beam  at  the  Laser  Window  Test  Facility  a window  of  uncoated  CVD 

zinc  selenide  was  exposed  for  10  seconds  to  a 2 kW  laser  beam  with  an 

approximate  gaussian  intensity  distribution,  A detailed  description  of 

(14) 

the  experimental  procedure  is  being  published  elsewhere.  In  brief, 

the  distortion  of  the  window  produced  by  heating  was  observed  by  taking 

50  frame  per  second  motion  pictures  of  Mach-Zehnder  interference  fringes 

during  and  after  exposure  to  the  laser  beam  using  a helium  neon  laser 

source  to  illuminate  the  interferometer.  During  and  subsequent  to  the 

exposure,  the  two  faces  of  the  window  were  surface  cooled  by  flat  jets 

of  nitrogen  gas  with  an  estimated  surface  heat  transfer  coefficient  of 
2 

0.  0124  W/cm  deg  C.  Other  pertinent  data  pertaining  to  this  experiment 
are  given  in  Table  I. 

By  analysis  of  the  motion  pictures,  a curve  of  fringe  count  versus 
exposure  time  was  obtained  fer  the  OPD  on  the  center  line  of  the  laser 
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TABLE  I 


EXPERIMENTAL  DATA  FOR  THE  CVD  Zinc  Se  TEST 
Window  diameter  2b  =3.6  cm 


Window  thickness 

Laser  beam  power 

Beam  mean  diameter 

Surface  heat  transfer  coefficient 

He-neon  laser  wavelength 

Nominal  absorption  coefficient 
( 12) 

Poisson's  ratio 
Elastic  modulus^ 

Index  of  refraction^ 

Temp,  coef,  of  n^^^^ 

Stress  birefringence^^ 

(13) 

Volumetric  specific  heat 

Thermal  diffusivity^^^^ 

Thermal  expansion  coef. 

(13) 


a 

. 82  cm 

p 

- 

2 kW 

l/a 

= 

. 24  cm 

h 

= 

.0124  W/cm^  deg 

A. 

= 

6.  33  X 10  ^ cm 

P 

= 

. 005  cm”  ^ 

V 

= 

.28 

E 

= 

6.  9 X 10^  ^ d/cm^ 

n 

s 

= 

2.403 

dn  /dT 

B 

= 

7.9x10  ^deg  ^ 

^11 

= 

-1.39  X 10  cm^/dyn 

*^12 

= 

0.  58  X 10  cm^/dyne 

pc 

= 

2.  65  J/cm^  deg 

K 

= 

. 049 1 cm  /sec 
„-6  -1 

a 

= 

7.2x10  deg 

a 

V 

= 

5.  6 X 10^  dyne/cm^ 

Yield  stress 


beam.  This  is  plotted  in  Figure  7.  Superimposed  on  this  curve  is  a plot 
of  the  theoretical  Eq.  (32)  evaluated  by  use  of  the  data  in  Table  I.  The 
window  was  supplied  to  us  by  AFCRl.  and  it  was  stated  that  the  previous 
measured  value  of  (3  was  0.  005/cm 

The  correction  term  to  the  t~ser  power,  for  surface  cooling, 
according  to  Eq.  (38)  is  evaluated  fir.«t.  After  9 seconds,  it  is  0.  985 
and  at  shorter  times  it  is  even  closer  to  1.0,  Therefore,  it  is  considered 
to  !)e  negligible. 

The  OPD  is  found  from  Eq.  (32).  The  first  term  in  the  bracket 

(thermal  expansion)  is  1.  18  x 10  ^ deg  \ The  second  term  (stress  optic 

effect)  is  13.  9 X 10  ^ deg  The  third  term  (change  in  index)  is 
-5  - 1 

5.  69  X 10  deg  and  is  the  dominant  one. 

Discussion.  Comparison  of  the  calculated  and  experimentally  mea- 
sured curves  of  Figure  7 shows  that  the  theory  accurately  predicts  the 
OPD  up  to  the  first  5 seconds  of  heating.  After  10  seconds  the  theututical 
curve  is  about  15%  below  the  observed  curve  and  has  a smaller  slope. 

There  are  three  possible  explanations  for  this  discrepancy.  The  central 
portion  of  the  window  has  been  heated  about  25  degrees  in  10  seconds. 

It  is  possible  that  some  of  the  properties  of  CVD  zinc  selenide,  besides 
n^,  have  a temperature  dependence  which  is  not  accounted  for  in  the  cal- 
culations. The  most  likely  property  to  have  a temperature  dependence 
is  P,  which  in  the  present  theory  is  assumed  to  be  constant.  Skolnik, 

ct  al.  , calculated  the  temperature  dependence  of  p on  the  assumption  that 

(15) 

it  is  due  to  multiphonon  absorption  and  stated  that  there  could  be  a 33% 
increase  in  P over  a 25  degree  temperature  change.  However  the  authors 
measurements  of  P for  the  samples  which  they  possessed  was  independent 
of  temperature  over  the  range  from  300  to  500  °C.  Therefore,  it  is  not 
likely  that  a P dependence  on  T is  the  correct  explanation. 

Another  possible  cause  of  the  discrepancy  at  long  heating  times 
is  that  the  present  theory  assumes  an  infinitely  large  window  with  a cen- 
trally heated  region.  In  actuality,  the  laser  beam  center  was  about  0.  7 cm 
from,  the  edge  of  the  window.  The  heated  zone  would  have  diffused  to  the 
edge  in  5 seconds  after  which  time  the  window  temperature  would  rise 
faster  in  this  experiment  than  for  a much  larger  window.  The  extent  of 
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Figure  7 Comparison  of  Experimental  and  Calculated  Optical 
Path  Difference 
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the  extra  heating  may  be  estimated  by  interpolation  of  Figure  1 at  ra  = 0 

2 

and  4k  a t = 22  and  44  corresponding  to  5 and  10  seconds  of  heating.  The 
temperature  of  a large  window  should  have  risen  by  50%  for  a 100%  increase 
in  the  time.  For  a smaller  window,  with  a smaller  diffusion  area,  after 
5 seconds  there  would  be  a smaller  increase  in  temperature.  This  is 
qualitatively  sufficient  to  explain  the  difference  between  the  observed  and 
theoretical  heating  curves  at  long  heating  times.  The  third  explanation 
is  that  the  value  of  a may  not  have  been  accurately  measured.  We  have 
calculated  that  for  the  conditions  of  the  experiment,  dT/T  = -3.  24  do /a  . 
Thus  a 15%  error  in  T could  be  accounted  for  by  a 5%  error  in  the  know- 
ledge of  n.  We  believe  that  a portion  of  the  discrepancy  at  long  heating 
times  may  be  ascribed  to  this  cause  also. 
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